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ABSTRACT 

Though pick-up ions (PUls) are a well known phenomenon in the inner heliosphere, their phase-space distribution nevertheless is a 
theoretically unsettled problem. Especially the question of how pick-up ions form their suprathermal tails, extending to far above their 
injection energies, still now is unsatistactorily answered. Though Fermi-2 velocity diffusion theories have revealed that such tails are 
populated, they nevertheless show that resulting population densities are much less than seen in observations showing power-laws 
with a velocity index of "-5". We first investigate here, whether or not observationally suggested power-laws can be the result of a 
quasi-equilibrium state between suprathermal ions and magnetohydrodynamic turbulences in energy exchange with eachother We 
demonstrate that such an equilibrium cannot be established, since it would require too high pick-up ion pressures enforcing a shock- 
free deceleration of the solar wind. We furthermore show that Fermi-2 type energy diffusion in the outer heliosphere is too inefficient 
to determine the shape of the distribution function there. As we can show, however, power-laws beyond the injection threshold can 
be established, if the injection takes place at higher energies of the order of 100 keV. As we demonstrate here, such an injection is 
connected with modulated anomalous cosmic ray (ACR) particles at the lower end of their spectrum when they again start being 
convected outwards with the solar wind. Therefore, we refer to these particles as ACR-PUls. In our quantitative calculation of the 
pick-up ion spectrum resulting under such conditions we in fact find again power-laws, however with a velocity power index of 
"-4" and fairly distance-independent spectral intensities. As it seems these facts are observationally well supported by VOYAGER 
measurements in the lowest energy channels. 
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1. Introduction 

Suprathermal ions, picked-up by the supersonic solar wind flow 
as ionized neutral atoms, have become known as pick-up ions 
(PUIs) and are produced all over the inner heliosphere with a 
typical upwind-downwind asymmetry wit h respect to the inflow 
direction of the neutral ISM inflow vector dRucinski et al.llT993t 
lFahr&Rucinskill999l) . In the case of PUI protons, their produc- 
tion is due to ph otoionization and charge exchange of interstel- 
lar H-ato ms (see Rucinski & Fahr 1991; Fahr & Rucinski 1999"; 
iRucinski et al. 2003; Bzowski et al, 2007 ). Their spatial distri- 
bution seems well understood, while the PUI phase-space trans- 
port is a much less settled subject. Especially it exists an ongo- 
ing debate of how efficiently pick-up ions just after the pick-up 
process are accelerated to higher energies due t o nonlinear wave- 
particle interactions (se e e.g.iIsenberg^l987;,Bogdan et alJl991 ; 
Fichtner et al. 1996; Fi chtneij|200 U IChalov & Fahjl 1 996l1 1 99~ : 



tion of these newly 


produced PUIs* is toroidal and unstable 


(see Winske & Lero\^ 


19841; IWinske et alJll985HLee & IdII 19871; 


Fahr & Ziemkiewicz 


19881). With the free energy of this un- 



Chalov et al.ii2004 ) and whether at all energy diffusion plays a 



relevant role in this transport. 

Some hint is given by the solar wind proton temperature be- 
havior with distance. The observed non-adiabatic temperature 
behavior namely proves that a specific solar wind proton heating 
must operate in the outer heliosphere which can only be due to 
energy absorption from pick-up ion generated turbulence, since 
con vected turbulence amplitudes quickly die out with distance 
(see lSmith etlDIMlllFahr & Chasheill2002i . 

Freshly injected PUIs represent keV-energetic protons in 
the supersonic solar wind frame and may be called here: 
"primary pick-up ions" (or: PUIs*). The velocity distribu- 



stable distribution PUIs* drive Alfvenic wave power. The lat- 
ter enforces p itchangle isotropization of th e initial velocity dis- 
tribution (see lChalov & Fahiill998l Il999al) . Due to wave-wave 
coupling, the wave energy generated by PUIs* at the injection 
wavelength A, = Us/Q.p is diffusively transported in wavevec- 
tor space both to smaller wavelengths where it can be absorbed 
by solar wind protons and to larger wavelengths where it is re- 
absorbed by all PUIs. This effect is seen as the main reason 
of solar wind proton heating occuring in the outer heliosphere 
(Smith et al. '2001; iFahr & Chash ei' 2002b; 'Ch ashei et al.]|2003l; 
Staw icki 2004). Only a small fraction of about 5 percent of the 
PUI-generated wave energy reappears in the observed proton 
temperatures. Freshly injected PUIs excite turbulences that can 
organize a power-law distribution. From this distribution, both 
the solar wind io ns and the PUIs the mselves can absorb en- 
ergy as shown by jChashei et al.l (l2003i) . Also the approach by 
llsenberg et alj ilQQ^ where energy diffusion of pick-up ions is 
not taken into account shows that only a low degree (2-5 per- 
cent) of the pick-up ion driven wave energy is absorbed by solar 
wind protons in form of thermal energy. This raises the question 
where the major portion of the wave energy produced during the 
primary pick-up process goes to. To clarify the energy redistribu- 
tions, kinetic and spectral details of the relevant processes have 
to be investigated. A detailed numerical study of the PUI ve- 
locity distribution and the spectral Alfvenic/Magnet osonic wave 
power evolution has meanwhile been carried out (Chalov et al.l 
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|2004 l2006bh and presents a simultaneous solution of a coupled 
system of equations consistently describing the isotropic veloc- 
ity distribution function of PUIs and the spectral wave power 
intensity. 

As one can see from this study, the largest portion of the 
self-generated wave energy is reabsorbed by PUIs themselves as 
a result of the cyclotron resonant interaction and leads to PUI- 
acceleration. It could perhaps be hoped that this energization of 
pick-up protons due to Fermi-2 stochastic acceleration processes 
eventu ally leads to the ubiquitous pow er-law PUI-tails pointed 
out by iFisk & Gloeckleil (l2006l l2007h . To the opposite how- 
ever, as reflected in the results presented bv lChalov et alj (|2004 
l2006b.) it is evident that this is not the case: even at larger dis- 
tances close to the termination shock (100 AU) the PUI distri- 
butions show a rapid cut-off" at energies higher than the injection 
energy. The question thus is raised here why power-laws have 
been seen at all. An explanation that we are favoring here is a 
new injection source to the PUI regime from high energies con- 
nected with modulated anomalous cosmic ray particles. These 
protons are primary ACR particles that occur with a spectrum 
down to the typical energy of the usually assumed PUIs. At 
this part of the spectrum, both particle species cannot be dis- 
tinguished. Therefore, we refer to them as ACR-PUIs. 

In Sect. |2] we investigate the physical possibility of power- 
law ions in the outer heliosphere as they are recently proposed 
by several authors and we find that they cannot occur with a 
power-index of -5. As we show in Sect. [3] the proposed pro- 
cesses are not effective enough to produce the desired ion tails 
which means that another mechanism has to lead to the observed 
spectrum. In Sect. H] we show how a high-energy source can be 
derived by taking a modulated ACR spectrum upstream of the 
solar wind termination shock. This injection mechanism is dis- 
cussed in Sect. |5] where we show that these high-energy ions 
can lead to power-law ion tails, however with a power-index of 
-4. The results are discussed and compared with observations in 
Sect.|6] 

2. Can power-law ion distributions be in equilibrium 
with hydromagnetic turbulence? 

Challenged by recent results concerning ion spectra at large dis- 
tances (McDonald et al. 2003; Decker et al. 2005; Kiraly 2005; 
iFisk & Gloeck ler 2006), we look into the problem of PUI phase- 
space transport under thes e new given auspices. Fh'st we dis- 
cuss the argument given by iFisk & Gloeckleii (l2006l l2007h that 
pick-up ions under resonant interaction with ambient compres- 
sive turbulences enter a quasi-equilibrium state with saturated 
power-law distributions of a somehow sacrosanct spectral ve- 
locity index of y,, = -5. 

The related, well known Kolmogorov formalism is based on 
a "dimensional" reasoning: The problem concerning the energy 
distribution in eddies of a typical wave number k is considered 
using two different dimensional quantities: namely the spectral 
energy density and the wave number k . The spectral energy 
flux is then defined by 



S = {kEk)/n 



(1) 



and has to be constant for stationary cases, is the typical life 
time of eddies with scale = 2n/k and is given by 



Tk = (kVk) 

with the typical vortex velocity v^^ given by 

vk = (Ekky^ 



(2) 



(3) 



Requiring a constant energy flux S given by Eq.([TJ, with 
Tk given by Eq.© then automatically leads to the Kolmogorov 
spectrum in the form Ek ~ k^^^^. 

In their thermodynami c al approach to particle spe ctra in 
equilibrium with waves, ,Fisk & Gloeckleii (l2006i l2007h . how- 
ever, only consider one single dimensional quantity, i.e. the ki- 
netic energy T, because in their case the energy distribution F 
has the dimension [F] - [T] ^ The flux combination 



AT 

Sf^TF — 
At 



(4) 



is then expected by them to be constant and, thus, the energy 
gain is given by 



(5) 



Ar 

At 

and hence the above requirements result in the requirement that 

F ~ 

This shows that the analogy of the Fisk-Gloeckler approach 
with the Kolmogorov formalism is not complete, since in their 
theory the quantity At is not specified. Only if we play the same 
game with the Kolmogorov turbulence using instead of Eq.(|2| 
the assumption t = (kuy^, where u is some independent external 
speed (especially independent on k), we will find from Eq.([T]) 
also the result 



Ek ~ k 



-2 



(6) 



analogous to the result obtained by Fisk-Gloeckler 

The problem may be briefly inspected here whether or not 
some spatial/temporal disturbances in the solar wind plasma can 
be considered as waves. They should be considered as waves, 
if the convection time Tconv = r/U is much greater than the 
passage period of these waves in the wind reference frame 
= 2jT/(k()VA.) where r denotes the radial distance to the sun 
and U the solar wind bulk velocity. The Alfven speed is desig- 
nated as va with the corresponding wavevector ko for the tur- 
bulence. For an inequality of the kind <s: fconv, this then 
leads to the definition of the principal wave turbulence corre- 
lation scale given by < rvp^jU ~ 0.1 AU. A more ex- 
act definitio n of the turbulence correlation scale is given in 
the paper of Chashei et al.l (l2003j). This val ue is also in accor- 
dance with the value used bv iFahii (|2007|) for his estimation 
of the upper possible velocity border. Some measurements in- 
dicate smaller values for the turbulence correlation length tha n 
calculated by us in this paper (see e.g. iMatthaeus et al]|2005h . 
However, the following conclusions in our paper drawn on the 
effects of turbulent heating of suprathermal protons become even 
less promising for an MHD-equilibrium to be established, if 
smaller correlation lengths prevail. One can expect that at larger 
distances, r > 1 AU, the value for increases proportionally 
to r, since a similar depend ence is also vah d for the outer scale 
of turbulence A:,^' (see Chas hei et al.ll2003b . Here we estimate 
the maximum energy for protons resonating with these largest 
scales Ln, at distances of about 100 AU and find with Vmax - 
Q(100AU)-L„(100AU)/r(y„,x) - ^(1 AU)-L™(1 AU)/r(v„ax), 
where v^^^ , yiv^^^) and Q are the maximal ion speed, the as- 
sociated Lorentz factor and the ion gyrofrequency, respectively. 
Protons with these speeds have energies of about fimax ~ 1 GeV 
independent on distance. Thus, ions with E < /Smax can be ex- 
pected to be scattered by waves, whereas ions with E > /Smax 
can only be scattered by large scale velocity structures in the 
solar wind. 

Since in any case, however, power-laws seem to be a fact 
well supported by observations within a certain range of solar 
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distances and ion energies, we shall nevertheless take this finding 
as serious here as deserved and determine in the following the 
absolute spectral intensity of this PUI power-law distribution and 
its consequences. 

If the PUI-distribution is given in the form 



(7) 



fpniir, V) = /pui.O ■ — 



with the velocity power index - -5, then the PUI density is 
given by 



«pui(?-) = 4;r/p 



pui,0 



VQ 



v^dv ^ 2;r/pui.o ■ 



(8) 



where /pui.o = /pui,o('") is a local normalization value, and vo and 
Voo are lower and upper velocity limits of the quasistationary PUI 
power law. 

Using t// = v/vo, one obtains for the power-law PUI pressure 
as the second moment of the distribution function the following 
result 



-Ppui('") = 4;r/pui,o 



'■^vlj iA'-V^#' = 27rm-/pui,ov^ln(iAoo) (9) 



Evidently, the definition of Ppui('") requires the determination 
of all three local values /pui,o('")> vo(r), and Voc{r) which we aim 
at below. 



2.1. Determination of tlie inner and outer velocity border 
Using the definition of the lower velocity border as done by lFahj 



(10) 



(120071) one finds 

1 1/4 



VO 



1 / r 



With the above result one obtains the PUI distribution in the 
form 



pui,0 



1 



ripuiir) ■ 



-3/4 



(11) 



The local PUI density can be derived from the interplanetary 
H-atom density tiuir, 0), which depends on the radial distance 
r and the inclination angle with respect to the upwind axis, 
and the eff'ective (charge exchange + photoionization)- induced 
injection rate ySpui = nii(r,0)[nsir)o-^,iU + Vi]. Here tisir), o-g^, 
U, v\ denote the solar wind proton density, the charge exchange 
cross section, the solar wind bulk velocity and the photoioniza- 
tion frequency. The H-atom density at distances r > 5 AU in the 
upwind hemisp here is satis factorily well given by the following 
expression (see lF^fTOTll) 



riuir, 0) - nn.oo exp 



/?pui,0'"o^ 



Ur sin 6 



(12) 



and thus leads to the following PUI density (also see 
iFahr&Rucinskill 19991) 

r 

npui(r, 0) = (y) «pui.o + j Pv-^r\ 0y^Ar' (13) 



From the above one derives the following radial space 
derivative which lateron in this paper will be needed 

g«pui(r, 0) _ -2npui(r, 0) /3j,ui{r,0) 
d? ~ ~r U 



(14) 



To calculate the PUI pressure, one more quantity in addi- 
tion is needed, namely the upper velocity bo rder Vco. To de- 
termine Voo w e follow the idea presented by Fisk & Gloecklej 
(I2006ll2007h assuming that PUI power law distributions result 
from a specific quasi-equilibrium state self-establishing such 
that the wave field transfers per unit of time as much energy 
to puis by energy diffusion, as energy is expended in the solar 
wind frame for the work done by the pressure gradient of co- 
moving puis against the magnetosonic fluctuations. The impor- 
tant restriction to energy diffusion by nonlinear interaction with 
the compressive fluctuations is that the typical diffusion period 
Tdiif ^ 3L^/v/i|| should be much larger than the con vection pe- 
riod given by Tconv - L^/U (see lChalov et al.ll2003l) . This leads 
to the requirement 



lOii 
3U 



(15) 



with as the mean free path for particles parallel to the mag- 
netic field. The uppermost velocity Voo is the limit at which this 
condition is just fulfilled: 



3UL^ 

An 



(16) 



We assume the interaction of the particles with a slab Alfenic 
turbulence field. The mean free path is given by 



^3v ni 

8 j D.' 



dju 



with the pitch-angle diffusion coefficient 



D 



1111 



(17) 



(18) 



from IChalov et alj (|2003|) for cyclotron resonant wave-particle 
interaction with unpolarized, one-dimensional, and isotropic tur- 
bulence, which leads to the velocity independent expression 



3/4 



for the mean free path, with the reference value 



D,,,o = 



E rE 



U 9L„ 



(19) 



(20) 



for the diffusion coefficient. Therefore, the upper velocity border 
is given by 



3UL^ 

Vco - — , 

which evaluates to 

9 /rE\3/4 



~ 0.3 I r / 



U 



(21) 



(22) 



by taking = 3 AU and ^||,e = 0.3 AU dChalov & Fahdl999bl) . 
The ratio 



— = 150x^' 

Vo 



(23) 
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is needed in the later calcul a tion a nd should be compared with 
the result obtained by iFahri (l2007h deriving the upper velocity 
border from the study of the upper-most resonance possibilities 
of ions with the largest prevailing correlation lengths existing 
in the solar wind velocity structures, yielding the result 



i/^oo = 30 ■ X 



-3/4 



(24) 



which gives smaller values than those derived for condi- 
tions when balanced^ pressure equilibrium in the sense of 
iFisk & Gloeckleii (l2007b is adopted. Here below we shall 
demonstrate that this value for ^oo, which is the direct conse- 
quence of the assumptions by Fisk, definitely leads to unrea- 
sonable consequences when we investigate the associated PUI 
pressure. 

2.2. The upstream pick-up ion pressure 

We calculate the PUI pressure resulting from power-law dis- 
tributed puis upstream of the shock and find with Eqs. (|9]l, ( fTOl i. 
and ^ 



(25) 



where we have introduced x - r/rE- We obtain the pressure gra- 
dient by differentiation: 



dP, 



dx 
dn 



= mU X 



dx ^ '2 

which leads with Eq. ( fT4l i to 



dx 



- mU X 



(26) 



(27) 



-^npui^"'^^ In (150jc"') + ^x^'^ In (l50x"') - «puiX 



.-1/2 



If now we derive the effective upstream Mach number, ne- 
glecting thereby solar wind electron and proton pressures com- 
pared to the PUI pressure, i.e assuming Pg; P^ « : Ppui, and fol- 
lowing the definitions by Fahr & Rucinskii (Il999l) . we obtain 



Mr = 



|(p, + p, + Pp„i) 

2ns 



(28) 



Inp^iX^/^ In (ISOx-^) - ^x3/2 In (150x-i) + n^^iX^/^ ' 

With ^ 10~^ and the PUI abundance f = n„Jn, ^ 0.2 

(see lFahr & Rucinskill 1999b . we obtain 

2X-1/2 



f ( I In ( 1 SOjc- 1 ) -I- 1 ) - A X In ( 1 50x- 1 ) ' 



(29) 



For a shock position at x ^ 100 and with the lower reference 
distance ro - 5 AU, this expression yields 



M*/ = 0.83, 



(30) 



which means that the effective upstream Mach number M* - 
0.91 is lower than 1. 

Hence, power-law PUI pressures in balance with the wave 
fields would not allow for the occurence of a termination shock 



of the solar wind (i.e., needing upstream Mach numbers M* > 
1!) 

The non-existence of the above "turbulence - particle" - equi- 
librium (TPE) can also be concluded along a different line of 
argumentations together with the calculation of the downstream 
PUI pressure connected with power- l aw dis tributed PUIs. Using 
the results presented in iFahr & Lavl 42000), one obtains on the 
basis of Liouville's theorem and conservation of magnetic mo- 
ment at the passage from upstream to downstream over the shock 
the PUI distribution function downstream of the shock (indicated 
by index 2) by the following relation 



fpui.li^) - —pfpm,l 



V 



(31) 



depending on the upstream distribution function /pui,i, where s 
denotes the compression ratio at the shock. From that relation 
one obtains for the downstream PUI pressure Ppm,2 the following 
expression 



pui,2 



2 



J /pui,i(v/ Vs)v'*dy 



Vsvo 



= ^« I /pui,l(v)vMv = S^Ppui,! 



(32) 



which states that the downstream PUI pressure is enhanced with 
respect to the upstream PUI pressure by the factor s^. Reminding 
that the latter pressure for TPE-conditions is given by 

nui.i('-sh) = »Jnpui(rsh)t/'(^j ln(l50^j (33) 



one finds that the downstream PUI pressure should amount to 

,1/2 

^pui 



^'pui,2('"sh) = OTnpui(rsh)i^f/^ I — 



rE 



In 1 150— 



(34) 



which normalized with the upstream solar wind kinetic energy 
density e^in - (l/2)mf/^(ns) would require that 



^'pui,2('-sh) ^ 2/'"sh 

= 2^(rsh)s — 



1/2 



In 1 150— 



(35) 



with ^(rsh) - 0.2 denoting the PUI abundance at the shock. This 
would mean that the downstream thermal energy of the PUIs 
is much higher than the kinetic energy of the upstream solar 
wind which is forbidden by physical reasons. This again leaves 
to conclude that PUIs cannot exist in pressure equilibrium with 
the compressional magnetosonic turbulence that was assumed 
by [pisk & Gloeckler (2007) and cannot be responsible for the 
stochastic particle acceleration up to regions near the solar wind 
termination shock. 

Requiring that the downstream thermal energy of the PUIs 
stays below the upstream kinetic energy of the PUIs would re- 
quire instead an upper border Vc» of the PUI power spectrum 
defined by 



pui,2 — 



J"/pui,i(v)v'*dv 



47!-/puiyVo ln((/roo) 



< yt/^(n.s + npui) 



(36) 
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yielding 

ln(l/r„) < 



U^(n^ + Hpui) f/^(ns + Wpui) 25 / r^h 



-1/2 

— I (37) 



Evaluating this formula with s - 2.5, ^pui = 0.2, and r^h - 
lOOrE leads to the result 



ln(i/roo) < 



25 



2 -2.52 ■0.2- 10 



1 



(38) 



and means that i/Zoa - exp(l) = 2.72 and Voo - 2.72t/, i.e. much 
lower than required for pressure equilibrium conditions. 

3. The relative effectiveness of energy diffusion and 
convective changes at ion phasespace transport 

In the following we want to clarify the role of energy diffusion in 
determining the shape of the PUI distribution function. We start 
from the transport equation adequate to describe the phasespace 
behavior of the PUIs by a distributio n function f(t, r, v) (see e.g. 
IIsenberg|[T987l:IChalov & F^ll996h 



+ Qir,v) + S(r,v) 



(39) 



where U is the solar wind bulk speed, Z),.,. is the velocity diffu- 
sion coefficient, and, Q(r, v) and S (r, v) are functions describing 
PUI- injection sources and - phasespace losses. Terms on the left 
hand side of the above Eq. ( |39] | under steady state conditions in- 
duce changes with a typical convection time Tcoiw given by 



r/U 



(40) 



where r is of the order of the heliocentric distance. 

Considering quasi-linear velocity diffusion of particles due 
to Fermi-2 type interactions with the Alfvenic- or magnetosonic 
turb ulence, both of which are leading to analogous expressions 
(see 



Toptvgin|[T98l iLe Roux & PtuskinI [1991 IChalov & Fahd 



For our estimates here, one can use the diffusion coeffi- 
cient derived by Schli ckeised ( Il989h which for estimate purposes 
can be represented in the following form 



^res 



o-l 



(41) 



where A^,es = Q./vis the resonant wave number, Q. the proton cy- 
clotron frequency, va the Alfven speed, ko the turbulence outer 
scale, a ^ 5/3 (or ^ 3/2) the power exponent of the ID- turbu- 
lence spectrum, and 6mo is the fractional turbulence level 



3 MO 



<B>2 



(42) 



where B is the induction of local interplanetary magnetic field. 
Typical velocity diffusion times Tdiff characterizing the action of 
the first term of the right-hand side of Eq. ( |39] | in changing the 
distribution function / can be defined by 



diff 



D,,v 



-2 



(43) 



Combining the relations given by Eqs.(l40]i. ( l42b . dTO . and ( l43T l 
yields as a typical ratio of the characteristic times 



X(v, r) = 



Tdiif 



3M0 



,,2 \ 



vU) 



(kresr) 



a-l 



(44) 



At heliocentric dista nces r ^ - I AU, one can assume 
dChashei et al.ll2003h 6mo ~ 0.1, va ~ O.IU, ko ^ lO ^cm"', 
or = 3/2 for Iroshnikov-Kraichnan turbulence (i.e., a power law 
for the spectral energy density ~ k^^^^) or a = 5/3 for 
Kolmogorov turbulence (i.e., Ei^ ~ k^^^^), which are the two 
mostly found and disc ussed spectral energy distributions in so- 
lar wind turbulence (IBale et alj2005l:lHorburv et alj2005h . Then 
near 1 AU, we find a ratio x(U, 1 AU) ^ 1 showing that con- 
vection and diffusion processes here are of comparable impor- 
tance for the particles with v ^ U, while convection effects 
are of increasing importance for particles with higher velocities, 
i.e. with V > U. This result is in good agreement with other con- 
siderations (iMatt haeus & Goldstein 1986; Jokipii & Kota 19891: 
Zank & Matthaeus 19921 IZanketal] 119961; iHorburv & Baloghl 
2001,) . ilsenbera (i2005i) . however, concludes that the effect of 
energy diffusion is small for PUIs. His similar treatment is 
parametrized by the percentage of turbulence occuring in slab 
Alfvenic fluctuations. Nevertheless, these two approaches can be 
conciliated if this percentage is assumed to amount about 20%. 

At larger heliocentric distances, r > 1 AU, the frozen-in 
Parker magnetic field near the ecliptic is nearly azimuthal and 
decreases with (1/r). Correspondingly, the resonant wave vector 
behaves like ^res ~ Q ~ {B } ~ which leads to ^res'" = const. 
Besides that, as argued by IChashei et alj (l2003h . the following 
r-dependences can be expected: ~ ' and Smo ~ r^^^^- 
Consequently, the radial dependence of ;^f(v, r) is exclusively de- 
termined by the r-dependence of ^mo meaning that 



-V/3 



(45) 



in case of Kolmogorov turbulence. The above Eq. ( |45] ) shows 
that at larger solar distances energy diffusion of particles cannot 
determine the shape of the suprathermal tail of the resulting PUI 
distribution function. 

In fact, assuming for reasons of a better clarification, just to 
the contrast here, a dominance of the energy diffusion, i.e. ne- 
glecting all the terms in Eq. ( |39] l except for the diffusion term, 
would deliver as a result of the transport equation the distri- 
bution / ~ V This distribution thus should evidently be 
much flatter than the ob viously observed distribution / ~ v"^ 
(iFisk & Gloecklerll2006 !) indicating that energy-diffusion in fact 
plays an inferior role. 

This result also turned out from more quantitive calcula- 
tions of wave-partic le i nteractions in the oute r heliosphere by 
IChalov et all (l2006al) or lFahr & Chasheil ( l2007h in which a con- 
sistent treatment of pick-up ion transport in phasespace and wave 
turbulence transport in fc-space was presented. From these calcu- 
lations it becomes evident that in fact in the outer heliosphere 
pick-up ions drive turbulence powers, but are not efficiently 
enough profitting from energy diffusion to produce power-law 
tails. 

The only way to fill up ion tails at energies larger than 1 keV 
resulting in some power-law distribution, as it seems to us, is 
to inject ions into the PUI regime at higher energies of about 
>50 keV and then let them cool to smaller energies at their co- 
convection with the solar wind. This would be a process similar 
to the one of pick-up ions originating from freshly i onized neu- 
tral atoms which are injected at about 1 keV (see iFahJ 120071 : 
ISiewert & Fahrll2b08h . In the following, we study such a high- 
energy injection mechanism to the PUI energy regime resulting 
from ACR ions that are adiabatically cooled to lower energies 
and, therefore, can be called ACR-PUIs. 
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4. Injection to the PUI regime from the high-energy 
side 

Here we first want to study the modulated part of the anomalous 
cosmic rays (ACRs) with origin near the solar wind termina- 
tion shock and investigate whether they can perhaps serve as a 
possible injection seed from high energies into the pick-up ion 
regime. We start from the analytic expression for the modulated 
ACR spectral intensity given for the case of a spherically sym- 
metric solar modulation. The differential ACR intensity for thi s 
case is given by the following expression dStawicki et al.ll2000h 

j(r,p) = p^f{r,p) 



1 s M-Hy-5 



U 

X exp ) 



(46) 



where S(rQ,pQ) denotes the ACR source function with source 
coordinates ro and po, /„ is the modified Bessel function of the 
first kind, U denotes the solar wind bulk velocity, and h and (p 
have the following definitions 



\'"0/ \PQ 
p \^ 



= 1 - 



\P0 



with the following additional quantities 

V roU 



yo 



(1 + 7 - (5)2 Krrfi 

2 + y 



(i+r-s) 



I 2 + y 

V - \l +y -S + — — e 



(47) 
(48) 

(49) 
(50) 
(51) 



and the following assumptions of the r- and p-dependences: 



U(r) ~ r^, Kyy{r,p) = /c„,o — — 

r^l \PE 



(52) 



with an arbitrary reference momentum /7r for the bulk velocity 
and the spatial diffusion coefficient a:,,.. 

As shown by Stawicki et al. ( 2000), the above expression can 
be simplified for low values of ACR particle momenta p and then 
yields the following expression for the ACR distribution 



f(r,p) 



bT{b) 



PoU 



(53) 



The modulation of the ACR spectrum due to spatial diffu- 
sion is contained in yp. This shows that, hence, the ansatz by 
IStawicki et al.l (l2000h is an analytical way to express the mod- 
ulated distribution function of the ACRs in the heliosphere. In 
the frame of its accuracy, this approach yields a valuable way to 
take the modulated ACR spectrum as the source of our cooling 
mechanism. The spectrum of the ACR particles at the solar wind 
termination shock serves as the ACR source function and, thus 



with g ood reasons, one adopts the function developed bv lDrurvl 
(Il983h for an ideally planar shock which is given by 

d 



1 . , 
S(ro,po) = -r—rQ^o 



(54) 



5acr,o- — exp 



Po 



5{ro - Tsh) 



where the spectral index q is connected with the shock compres- 
sion ratio i by ^ = lisjis - 1). p\, p^^x are appropriate values of 
the Fermi- 1 injection momentum and of the upper cut-off mo- 
mentum, respectively. This then leads to 



f(Ap) 



3S 



CO 

ACR.0 r dpolpo 

bFibW j Pi [p, 



-9-1 



exp 



PO 
7^ max 



(1 + y-6y Krr,>ih 



exp 



(1 + 7 - 6y 



,sh 



(55) 



This expressions repres ents a constant distribution function. 
Numerical simulations bv Fichtner & Sreenivasan (1999) show a 
linear dependence of the ACR flux on the energy (Je ~ E). This 
means that the distribution function in phase-space is constant 
which can be shown by considerations about normalization of 
the differential flux and the distribution function: 



/ = 47r 



m jE dE 
P^ <ip 



jE 

E 



— - const. 



(56) 



Therefore, the f{r, p) of the ACR particles is given by a con- 
stant value which follows from modulation theory: 



f{r,p)^U 



(57) 



From the above expression we deve l op the total s tream- 
ing 5(r,pp„i.i) (see Gleeson & AxfordI Il967t iFahil 1 19901: 
iFahr& VerscharenI 120091) for the ACR-PUIs in flie range be- 
tween the injection momentum pi and an upper border p-pun- 
ACR particles with higher momenta propagate only according 
to the transport equation for cosmic rays and do not participate 
in the diffusive injection. With these considerations we find 



^ ] '%^~ypTpfr^ — 9 — 



(58) 



as the total particle streaming at each place r. Since we expect 
that Ppuii s> Pi, the lower border pi can be neglected in Eq. (|58] |. 

Now we consider the local spatial divergence of this ACR- 
PUI streaming at some lower momentum border and take this to 
be the high-energy injection source Q to the PUI regime. With 
this choice, we obtain 



^ 1 5 2/7 iifZ/cons, 



(59) 



as the divergence of the ACR particle flux, constituting a PUI 
source. 

5. PUI Boltzmann-Vlasov equation with 
ACR-induced injection 

In the "solar" rest frame (SF), the representative Boltzmann- 
Vlasov equation (EVE) is given by 



, dv 



V., / 



df dv 
+ U — 
or or 



df 
dv 



(60) 
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under the prevailing conditions in the outer hehosp here, i . e. neg - 
ligible energy diffusion, for the stationary case (l Fah3l2007l) . 
The corresponding BVE equat ion in the "solar wind" rest frame 
(WF) has the form (lFahiil2007l) 



df 1 d / J -X - 



(61) 



where the second term on the left hand side describes the 
velocity-space divergence of the phasespace flow connected with 
the magnetically induced deceleration (i.e. magnetic cooling), 
which is indicated by the subscript m. The coordinate t = f(r) 
denotes the proper time in the co-moving reference frame (WF). 
Furthermore, P{r, v) is the local ion injection rate to the PUI 
regime. 

For low energies (< 1 keV), this rate is due to locally 
freshly ionized neutral H-atoms and is given by P{r, v) - 
P(r):^^pjj S(v - U) with I3(r) b eing the local pick-up ion produc- 
tion rate jFahr & SiewertlBOOSI) . 

For high energies (> 50 keV for the ACR-PUIs), the relevant 
injection, in contrast to the normal sub-keV PUIs, in our view is 
given by the upper expression derived from the modulated ACR 
spectrum and given by Eq. i59[ . 

Furthermore, v'm is determined by the magnetically induced 
velocity decrease of particles with a velocity v, when they are 
convected outwards with the solar wind bulk flow at a mean ve- 
locity U to larger distances where the co-convected interplane- 
tary magnetic field B appears reduced in magnitude (Fahr 2007). 
At larger distances r > 5 AU near the ecliptic, the magnetic field 
decreases like (1/r) (i.e. in case of the nearly azimuthal, distant 
Parker field). Under these conditions, one finds the following v- 
dependent magnetic velocity-space drift 



V 

-U- 

r 



and its associated radial gradient 



dr 



1 

— i 

U 



(62) 



(63) 



in the form given bv lFahr & Siewertl(l2008h . 

Ions which are picked up at with a velocity U will, with- 
out other processes being involved, have "magnetically" cooled 
down to a velocity v at r if the relation r,.(v) = rvjlJ is fulfilled. 
The injection of freshly created pick-up ions at r,,(v) with an ini- 
tial velocity v - U will be responsible for ions with velocity v at 
r. 

Taking all these constraints together, one finally finds, when 
reminding that the time and distance coordinates are related to 
eachother by dr = f/df, that for velocities v <U the solution f or 
/ in the WF is given by (see lFayi2007t ISiewert & Fah^l 2008h 



1 



27r U 



This distribution function /< can be evaluated for larger so- 
lar distances r > = 5 AU in the upwind hemisphere for near- 
ecliptic positions assuming that at such solar distances the up- 
wind H-atom density can be considered as essentially constant, 
meaning that hh (77'") = n^iir) - hh.co . This in fact is an accept- 
able approximation for solar distances r > 5 AU and velocities 
lU >v>Q.2U, and then leads to 



/< 



2nU 



Vex.E'-E 



2:7rr 



i.e. to the astonishing fact that under pure magnetic cooling, the 
resulting PUI distribution function for velocities v < If/ is a 
power law with the interesting power index a - -5 predicted 
and confirmed bv iFisk & Gloecklej (l2006l l2007l) : however, in 
their case expected as result of an assumed quasi-equilibrium 
state established between magnetoacoustically driven ion energy 
diffusion and magnetoacoustic turbulence generation. 

The above result is valid only for ions with v <U . One possi- 
bility that we now start to see here is that for ions with velocities 
V > If/, i.e. much higher than the original PUI injection thresh- 
old, one has to consider in addition some high-energy injection 
rate due to modulated ACRs as we have derived above. 

For those ACR-induced ions (i.e. for the ACR-PUIs) with 
the relevant source function 



Q 



9? ' 



(66) 



we find analogously to Eq. ( l64l i for ACR-PUIs 



/> 



1 



9r 



2n 



U 



2n 



U 



(67) 



which finally leads to an r-independent distribution function for 
the range v > U in the form 



I 1 3 4 r -4 

/> = ^'"pVpui,i/co„stV 



(68) 



Alltogether, thus, we obtain the total PUI distribution in the 
following form: 



/ = /<+/> 

yex,ErlU 

2nr 



«H,ooV '•HilJ -v) + £-m^v;„i ./cc 



1 

9^ 



(69) 



where H(x) is the well-known step function with H{x < 0) = 
Oandi/(.x>0) = 1. 

In any case, the above result shows that a velocity power- 
index of "-5" is obtained for the velocity-range v < JJ < 
and a power index of "-4" is obtained for the velocity-range 
U <v < Vpui.i where Vpui,i = p^m,\lm^. 

Rescaling velocity in units of U (i.e. X - v/U), leads to the 
total distribution function 



(o-exnsEU) r?- 



1 



—-^nn,^X-'mi-X) + -,nX^^J,,„,,X-' 



2nU^ 



A 



(7) 



X-^H(l-X) + 



2 mXuuU'U 



9n 



ns,E 



-X-"* (70) 



We introduce now two dimensionless characterizing quanti- 
ties. First, we define 



(64) ^ = (o"exnH,<x,rE) - 2.3 x 10 



-3 



where the c harge exchange cross section o-g 
(Fahr 197lh and the hvdrogen density mh.oo 
(llzmodenov et al.ll2003l) are used. 
With 



(71) 
0.15 cm-3 



10-1^ cm^ 



A = 



2 'nlxU^U'f.. 



(65) f(r,v) 



97T «s,E 

one can finally find 

nsE 



2kU^ 



A 



X'^H(1-X) + A-X'* 



(72) 



(73) 
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ACR-PUIs 
PUIS and ACR-PUIs at 70 AU 




Fig. 1. Phase space density spectra for PUIs and ACR-PUIs. The 
PUI distribution is dependent on the distance from sun, whereas 
the ACR-PUIs are not. PUIs cannot gain velocities above U, 
i.e. X - \. The solid line is led down to the ACR-PUI spectrum 
at the cut-off at X = 1 artificially. 

Measurements by the SWICS instrument on the ULYSSES 
space probe show a PUI phase space density of about 
lOOs-'km"^ at two times the solar wind velocity at a distance 
of 5.26 AU (lGloeckleill2003l) . Eq. ^ leads to a value for the 
PUI phase space density of about /< ^ 155 s^ km"'', which is in 
good accordance to the observations. The difference could be a 
consequence of the assumption of a constant n-n, especially at 
small distances from the sun. 

Observations by the L ECP instruments on VOYAGER 1 and 
2 jLanzerotti et aLlKOO l*) show an r-independent proton inten- 
sity in the energy range 0.57-1.78 MeV of about - A-x 10""^ 
particles per (cm^ s sr MeV), which corresponds to Je - 250 
particles per (cm^ s sr erg) in Gaussian units. This value can be 
converted to a phase space density (cf. Eq. (|56]l): 



/> =4;r 



p\\ MeV) 



2.7 X 10"^'' s^ cm-*' 



(74) 



If we assume that these ions observed bv' Lanzerotti et al] (1200 lb 
are ACR-PUI particles with the proposed behaviour, this leads 
to a value for A obtained from the measurements of 

A = ^^^4(1 MeV)/> = 3.1 X 10"'°. (75) 

We take this observational value to avoid uncertainties at the de- 
termination of the ACR intensity /const and the injection border 
-'^pui.i- The calculated spectra are shown in Fig.[T] 

6. Comparison of our results with data 

Now we want to compare the above results with Voyager and 
Ulysses spectral ion data. For that purpose, we first transform the 
above velocity distribution function into a distribution of spectral 
energy flux j{E) given in units of [ions/(cm- s sr MeV)] . This 
then leads to 

\(l-»,0/2 

j{E) ~ VI-/(£) 



[E'oI 



(2-a,)/2 



where a,, and oe are the velocity- and the energy -power indices, 
respectively. As one can easily see, av - 5 leads to o-^ = 1.5, 
whereas a,, = 4 would lead to oe = 1.0. Looking into the data 
obtain ed by VOYAGER 1 during the year 2004 (see lDecker et al.1 
120051 12OO 6) before the shock crossing occured, one can find - 
even though the data are very time-variable in this intervall - that 
for ions with higher energy (E ^ 10^ keV) the energy power in- 
dex has been observed with aE(E > 200 keV) ^ 1.5, whereas in- 
terestingly enough for ions with lower energy (i.e. £ ^ 40 keV) a 
smaller power index of aE{E < 200 keV) ^1.0 seems to be in- 
dicated. That supports our above derived theoretical prediction 
that the velocity power index of ACR-induced pick-up ions in 
the outer heliosphere is o-,. ^ 4, rather than a,, ^ 5 as for the ions 
with higher energy (MeV). 

Also duri ng the mo st recent VOYAGER-2 crossing of the 
shock (see De cker et al .' 2008) it became evident that the energy- 
power indices registered both before and after the shock crossing 
show values of aE,i - aE.i - 1.2 ± 0.2 which also nicely con- 
firms a velocity index of a,, ^ 4. 

In addition, it is important to recognize that the particle de- 
tectors of the VOYAGER spacecraft did not see radial changes 
of the spectral flux intensity jiE) at the lower energy channels 
0.35-1.5 MeV) during the period frorn years 1995 through 2005 



iLanzerotti et al teOOlMKrimigis et al.ll2003l:lDecker & KrimigisI 
I2OO3I) . This is also a support for the newly derived theoreti- 
cal expression given in Eq. ( l73l l and showing that a distance- 
independent intensity can be expected at energies 5 keV < E < 
100 keV. Only the part of the spectrum that is a result of cooled 
pick-up ions shows a dependence on the distance due to the dif- 
ferent PUI abundances. The ACR-PUIs that are injected from 
ACRs and cooled from high energies, however, are proven to 
have no radial variation. 

Our model explains occuring power-law tails with an index 
-4. There are, however, observations around 5 AU undoubt- 
edly show ing a power-index -5 for suprathermal ion tails (e.g., 
iGloeckler 2003^ As we can show, the favored explanation by 
Fisk & Gloec k!erl(l2006ll2007l) cannot hold in the described way. 
Maybe other processes (such as anomalous Fermi-2 type wave- 
particle energy diffusion), which are not treated in our approach, 
can lead to a power-index of -5 at energies just above 5 keV. 
But this interaction does not expla in the occurence of extended 
power-law tails at larger distances jChalov et al.ll2004l) . 

It is perhaps still a little bit an open question, where, i.e. at 
what lower energy, one should cut off the modulated ACR spec- 
trum to calculate the injection to the PUI regime. Even though 
this does not count very much in quantitative terms, since the 
lower energy part of the modulated ACR spectrum has in fact a 
constant distribution function (see Eq. (|56]l), it may nevertheless 
represent an intellectually interesting question, to decide up to 
what energies convection and Fermi-2 energy-diffusion are dom- 
inant, and from what energies upwards Fermi- 1 acceleration and 
spatial diffusion dominate. There is one clear hint given to an- 
swer this question, namely by the momentum-dependence of the 
spatial diffusion coefficient Krr{r, p) generally given in the form 



K(r, p) = KoKr(r)Kp(p) ~\—] { — 



(77) 



(76) 



where, dependent on the momentum-range , the exponent y is ex- 
pected to be in the range 1 < y < 2 (see I Le Roux & Potgietej 
Il992t |Jokipiil[T996l) . This shows that spatial diffusion becomes 
less and less efficient, the lower the ion momentum p is. Below 
some critical value p^^, ions loose their degree of kinetic free- 
dom to spatially diffuse relative to the solar wind background 
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flow and, thus, they are simply convected outwards with the so- 
lar wind bulk flow then. 

The determination of the absolute height of the ACR-PUI 
spectrum is a slightly problematic endeavor. In our calculation, 
we normalize the absolute spectrum with the aid of observational 
data. The ACR intensity /const and the limit for the momentum, 
up to which the injection is effective due to the total streaming 
(see Eq. dSSll), are quite uncertain. Another open question re- 
sults from a lack of data. Especially, the transition from the PUI 
to the ACR-PUI part has not been covered sufficiently by parti- 
cle detectors in the outer parts of the heliosphere. Maybe future 
missions can provide a closer look to this part of the spectrum. 
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